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In the present work we include the isovector-scalar 5-meson in the quark-meson coupling model 
(QMC) and study the properties of asymmetric nuclear within QMC without and with the 5-meson. 
Recent constraints set by isospin diffusion on the slope parameter of the nuclear symmetry energy 
at saturation density are used to adjust the model parameters. The thermodynamical spinodal 
surfaces are obtained and the instability region at subsaturation densities within QMC and QMC5 
models are compared with mean-field relativistic models. The distillation effect in the QMC model 
is discussed. 
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I. INTRODUCTION 

The instabilities presented by a system are directly re- 
lated with the possible phase transitions it can undertake. 
At subsaturation densities a liquid-gas phase transition 
in nuclear matter is predicted and it is normally tested 
in nuclear reactions. The formation of highly excited 
composed nuclei in equilibrium with a gas of evaporated 
particles can be interpreted in the framework of hydro- 
dynamics as two coexisting phases of nuclear matter, a 
liquid and a gas phase. During these reactions, phase 
transitions may occur depending on the temperature and 
densities involved. The liquid-gas phase transition also 
plays an important role in the description of the crust of 
compact star matter at densities between 0.03 frn and 
saturation density (~ 0.15 fm~ 3 ). It essentially consists 
of neutron rich nuclei immersed in a gas of neutrons. It 
has been shown that this phase transition leads to an 
isospin distillation phenomenon: the isospin content of 
each phase is different, most of the gas being composed 
of neutrons and the liquid being closer to symmetric mat- 
ter [3 

In the present paper, we employ the quark-meson cou- 
pling model (QMC) 0, Q to investigate the thermody- 
namical instabilities of asymmetric nuclear matter. In 
the QMC model, nuclear matter is described as a system 
of non-overlapping MIT bags which interact through the 
exchange of scalar and vector mean fields. An earlier 
study of ANM within this model has been focused on 
the effect of isospin asymmetry and temperature on the 
equation of state and on the coexistence surface [ij . We 
here consider an extension of the model that includes the 
scalar isovector virtual 5(a (980)) field ||. Its presence 
introduces in the isovector channel the structure of rel- 



ativistic interactions, where a balance between a scalar 
(attractive) and a vector (repulsive) potential exists. The 
5- and p-mesons give rise to the corresponding attractive 
and repulsive potentials in the isovector channel. The in- 
troduction of the 5-meson will affect both the behaviour 
of the system at low and high densities. In the last case 
due to Lorentz contraction, its contribution is reduced, 
leading to a harder equation of state (EOS) at densities 
larger than ~ 1.5 po [a]- At low densities a reduction of 
the symmetry energy will occur which will allow for more 
asymmetric matter. 

In 0, Q the instabilities in ANM have been investi- 
gated within relativistic mean field hadron models, both 
with constant and density-dependent couplings at zero 
and finite temperatures. It was shown that the main dif- 
ferences occur at large isospin asymmetry and at finite 
temperature. In particular it has been shown that the 
predicted density at the inner edge of the crust of a com- 
pact star, from the crossing of the /3-equilibrium equation 
of state (EOS), is model dependent Q. 

In the present work we investigate thermodynamical 
instabilities within the QMC model with and without 
the isovector-scalar 5-meson. Although in this model the 
isoscalar vector channel described by the p meson is in- 
cluded in a similar way to the non-linear Walecka model 
(NLWM), the non linearities in the a and 5-fields arise 
from the minimization of the bag energy. In particular, 
the NLWM used in [f| does not include non-linearities in 
the 5-meson. We may therefore expect a different behav- 
ior of asymmetric matter. 

The paper is organized as follows: in section II an ex- 
tension of the QMC model to include the 5-meson is dis- 
cussed, in section III we make a short review of the cal- 
culation of the spinodal surface, in section IV results are 
presented and discussed and some conclusions are drawn 
in the last section. 
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II. THE QUARK-MESON COUPLING MODEL The total energy density of the nuclear matter reads 



In what follows we present a review of the QMC model 
and its generalization to include the iso-vector-scalar 6- 
meson. 

In the QMC model, the nucleon in nuclear medium 
is assumed to be a static spherical MIT bag in which 
quarks interact with the scalar (a, 5) and vector (u>, p) 
fields, and those are treated as classical fields in the mean 
field approximation (MFA) [1, The quark field, ip qi , 
inside the bag then satisfies the equation of motion: 

[ip - H-gla-glrJ^-giuj 

+ 2 9 p TzPo31 ° ^n( x ) =0 ' <l = u,d (1) 



1 2 2 i 1 2 2,1 22 ,1 2c2 

e = ^/ + -m u w + -m p p 03 + -m s 6 3 



J2— / k 2 dk[k 2 + M* N 2 (a,S)] 



1/2 



N 

and the free energy density is given by 

J~ £ [LpPp PnPni 

where the chemical potentials are given by 



k 2 p + M; 2 + g^p- 



fJp 



P03, 



(8) 



where mi is the current quark mass, and g q 



9 q s, 



and 



9 

g q p denote the quark- meson coupling constants. The nor- 
malized ground state for a quark in the bag is given by 



i/) qi (r,t) = Af qi exp(-ie qi t/Ri) 
hi (x qi r/Ri) 



Xq 



4tt 



(2) 



Pn 



fc2 + M, 



" 2 + 9up- yPo3- 



The vector mean field ujq and po3 ar e determined 
through 



LOq 



9u>(p P +Pn) 9p(Pp-Pn) , a s 
5 > P03 — 5 ' v y J 



where 



£ g, = tt qi +Ri [ 9% u + ^9pTzPo3 ) ; Az« 



Q qi - Ri m* 



a, 



Ri m* 
(3) 



with the normalization factor given by 

M q 2 = 2Rfj$(x q ) [n q (n q - 1) + R N m* q /2] I x\ , (4) 



where Q„ 



:, ; =, ^x 2 qt + (R N m*) 2 , m* = m°-g$ a-g q 5 r z S 3 , 

Ri is the bag radius of nucleon i and \q is the quark 
spinor. The bag eigenvalue for nucleon i, x qi , is deter- 
mined by the boundary condition at the bag surface 



hi(x qi ) = PqijliiXqi) 



(5) 



The energy of a static bag describing nucleon i consisting 
of three quarks in ground state is expressed as 



E 



bag 



Ri 



Zi 
Ri 



itR?B n 



(6) 



where Zi is a parameter which accounts for zero-point 
motion of nucleon i and -Bat is the bag constant. The set 
of parameters used in the present work is given in Ref. 
hfj. The effective mass of a nucleon bag at rest is taken 
to be M* — E^ as . The equilibrium condition for the bag 
is obtained by minimizing the effective mass, M* with 
respect to the bag radius 



dR* 



0, i — p, n. 



(7) 



where g^ = 3g q and g p 
and 5 3 are fixed by 



9l- 



Finally, the mean fields <7q 



de 



0. 



de 



(10) 



In order to set the model parameters, we start by fix- 
ing the free space bag properties. They are obtained by 
fitting the nucleon mass and enforcing the stability con- 
dition for the bag in free space. We consider two sets of 
free space parameters, taking an equal proton and neu- 
tron mass value in a first moment, and then proceeding 
by considering different proton and neutron masses after. 

In the first case, we consider the bare nucleon mass 
M = 939 MeV and the bag radius, R p = R n = 0.6 



fm. The unknowns Z„ 



3.986991 and B 



1/4 
N 



211.30305 MeV are then obtained by setting the nucleon 
bag energies to that (single) bare nucleon mass value. 

In the next step, we take the physical nucleon mass 
values as M p = 938.272 MeV and M n = 939.56533 MeV 
and the bag radius for protons as R p = 0.6 fm. The un- 
knowns Z p = 3.98865, Z n = 3.98471, B 1 ^ = 211.26209 
MeV and the neutron radius R n — 0.6002 are then ob- 
tained. Note that for fixed proton bag radius R p = 0.6 
we observe a decrease on Z n and on the bag parameter 

1 /4 

B^ for the nucleons. Next, we fit the quark-meson cou- 



pling constants g%, g q s , g w = 3g q and g p = g* p 



and g p = g q p for the 
nucleons, so as to obtain the correct saturation proper- 
ties of nuclear matter, Em = e/p — M = —15.7 MeV at 
p = po = 0.15 fm~ 3 , a sym = 33.7 MeV. For the couplings, 
we have g q — 5.981, g u — 8.954. In our first case (when 
no effective mass difference between p, n is considered) , 
g p = 8.615. 
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TABLE I: Nuclear matter properties of the models used in the present work. All quantities are taken at saturation, except the 
density p s for which the pressure has a minimum and the incompressibility is zero. 



Model 


B/A 


Po 


K 


M* /M 


£sym 


L 






Ps 




(MeV) (mi" 3 


) (MeV) 




(MeV) (MeV) (MeV) (MeV) (fm" 3 ) 


NL3 [16] 


16.3 


0.148 


269 


0.60 


37.4 


118.3 


101 


-608.8 


0.096 


NL3<5 [18] 


16.3 


0.148 


270 


0.60 


37.4 


153.1 


427.1 


-491.5 


0.096 


TW [17] 


16.3 


0.153 


240 


0.56 


32.0 


55.3 


-125 


-456.8 


0.096 


QMC 


15.7 


0.150 


291 


0.77 


33.7 


93.5 


-10 


-570.8 


0.098 


QMC<5 


15.7 


0.150 


291 


0.77 


34.2 


102.1 


34.8 


-577.6 


0.098 


BHF [19] 


14.7 


0.182 


176.5 


0.79 


33.2 


63.4 


6.04 


-374.3 


0.119 



The properties of asymmetric nuclear matter have re- 
cently been related to both terrestrial data and star prop- 
erties from Vela pulsar glitches, which sets the symmetry 
energy slope value to L = 88 ± 25 MeV [H We 
then consider the <5-meson and determine the values of 
the couplings so as to have L = 102.077 MeV, which sets 
9pN = 12.599 and g q s — 12.6. In this case, the (p,n) mass 
splitting manifests in the different values for the effective 
masses: M* = 727.718 MeV and M* = 729.007 MeV, at 
saturation. 

We take the standard values for the meson masses, 
namely ra a — 550 MeV, m u — 783 MeV, and m p = 770 
MeV. 



in the latter, as it defines the spinodal surface, which is 
determined by the values of T, p, and y p for which the 
smallest eigenvalue of Tij becomes negative. The asso- 
ciated eigenvector defines the instability direction of the 
system, in isospin space. 

It has recently been argued [3] that in ANM the spin- 
odal instabilities cannot be separately classified as me- 
chanical or chemical instabilities. In fact, the two con- 
ditions that give rise to the instability of the system are 
coupled so that the instability appears as an admixture 
of nucleon density and concentration fluctuations. In the 
following we study the direction of instability and the 
spinodal for the different models considered. 



III. STABILITY CONDITIONS 



IV. RESULTS AND DISCUSSIONS 



The stability conditions for asymmetric nuclear mat- 
ter, keeping constant volume and temperature are ob- 
tained from the free energy density imposing that this 
function is a convex function of the densities p p and p„, 
i.e. the symmetric matrix with elements 



T = 



d 2 T 



dpidpj j T 
is positive 0, [H, O • 

This is equivalent to imposing 



(11) 



dp P ' d(p p ,p n ) 



(12) 



where we have used pi = 

The two eigenvalues of the stability matrix are given 

by m " 



*± = \ (Tr(^) ± VTr(-F) 2 - 4Dct(^)) , 
and the eigenvectors 6p± by 



(13) 



6pf 



J J 



I, j =p,n. 



The largest eigenvalue is always positive whereas the 
other can take on negative values. We are interested 



In the present section we compare the model properties 
of QMC and QMC5, respectively with and without the 
(5-meson, the non-linear Walecka model (NLWM) NL3 
[ill, with and without the 8- meson, and the density de- 
pendent relativistic hadron model TW [l7|. We will also 
refer to the nuclear matter properties obtained within a 
microscopic Brueckner-Hartree-Fock (BHF) approxima- 
tion using the realistic Argonne V18 nucleon-nucleon po- 
tential plus a three-body force of Urbana type [19( . 



A. Model properties 

We will first compare the equilibrium properties of nu- 
clear matter described by the different models consid- 
ered. The parameters of these models have been fitted 
to similar binding energy and saturation density values 
as seen in Table |T] At saturation, the effective mass in 
QMC is much larger than the corresponding mass in the 
other models, which is a characteristic of the model [![. 
In Fig. [TJ it is seen that the QMC mass decreases much 
slower with density. Even the hadronic models we study 
show quite different behaviour among themselves. NL3 
has an almost linear decrease on the mass whereas TW 
has much faster drop at low densities, and shows a less 
dramatic fall as density increases, crossing the curve for 
NL3 at p ~0.18 fm -3 . Incompressibility is one of the 
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FIG. 1: (Color online) Effective mass for QMC (brown, dot- 
dashed), NL3 (green, dashed) and TW (red) in symmetric 
matter. 



bulk properties that distinguishes the different models, 
but it is on the isovector channel that lies the largest dis- 
tinctions among the different models we use. Although 
having identical (or barely different) values for their bulk 
isoscalar properties, similar models differ considerably on 
the isovector parameters we discuss next. 

We now compare the symmetry energy and its slope 
and compressibility for all models in this work (Fig. [2] 
and Table [J) . The symmetry energy in our relativistic 
mean field models is given by 



6e F 2 



Am 2 



Mi 



(14) 



Mq 2 is the Fermi energy of the nucle- 



where ep = \JP 1 
ons, and Mq is their effective mass in symmetric matter. 
The NL3 models have the largest value at saturation of 
the models considered, 37.4 MeV. 

The symmetry energy slope L(p) is defined by L = 
ZpodEgym/ dp. The curvature parameter of the symme- 
try energy K sym = 9p^d 2 £ sym /dp 2 (Fig. Gfc)) is also 
of interest because it distinguishes between different 
parametrizations. In particular, the quantity K asy = 



6L can be directly extracted from measurements 



of the isotopic dep endence of the giant monopole reso- 
nance (GMR) |20|. Recent measurements of the GMR 
on even-A Sn isotopes give a quite stringent value of 
K aS y — —550 ± 100 MeV. According to this value, the 
hadronic and QMC models we use here (see Table Q| sat- 
isfy the above constraint, whereas the BHF results lie 
slightly below. 

The symmetry energy within QMC and QMC<5 shows 
an extremely linear behaviour with density (Fig. [2^)), 
in comparison with all hadron models shown. This is 
quite visible from the symmetry energy curves, but un- 
doubtedly clear from the slope parameter L (Fig. [2p). 
At larger densities the symmetry energy in QMC is es- 
sentially defined by the second term of Eq. I|14p. pro- 
portional to the density, due to the small variation of 
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FIG. 2: (Color online) Symmetry energy (a) and its slope pa- 
rameter L = 3po£' sym (b), and K sym (c) in the QMC5 (blue, 
dotted), QMC (brown, dot-dot-dashed), NL3 (dark green, 
dashed), NL35 (light green, dot-dashed) and TW (red) mod- 
els, for symmetric matter. 



the nucleon effective mass with density. Although still 
quite hard above saturation density, the QMC symmetry 
energy is softer than NL3, but harder than TW. At sub- 
saturation densities and considering only models without 
the (5-meson, the QMC symmetry energy takes the small- 
est values. The introduction of the 6 has the expected 
effect: at subsaturation densities the symmetry energy is 
softer but above saturation values it becomes harder due 
to the saturation of the (5-meson field Q . In Fig. [2J we 
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show both the NL35 and QMC<5 symmetry energies. The 
effect of the (5-meson on the QMC at subsaturation den- 
sities is quite small, much smaller than the effect seen 
in NL3. It is above the saturation density that the 8- 
meson has a larger effect in QMC. From the slope of the 
symmetry energy it is seen that while for QMC the slope 
decreases slightly with density, for the QMC<5 model it 
increases slightly, with a value close to 100 MeV. In the 
bottom figure we also plot K sym . The 5 has a very strong 
effect in the NL3 model. The QMC model is less affected 
but in both cases the presence of the 5 increases the sym- 
metry incompressibility K sym , becoming slightly positive 
for QMCfS while it was slightly negative for QMC. The 
model TW is presenting the smallest values. 

Fig. [3] shows the proton and neutron radii for proton 
fractions y p =0.5 and 0.0. For symmetric matter, the neu- 
tron bag is larger in this model due to the mass proton- 
neutron difference. This result has also been reported 
in 0. Decreasing the proton fraction increases the pro- 
ton radius and the neutron and proton radii cross at a 
certain value of the density, isospin dependent: density: 
-0.175 fm- 3 for y p =0.3, -0.11 fm~ 3 for y p =Q.l and -0.1 
fm~ 3 for neutron matter. The bag radius is sensitive to 
isospin content by a small amount. It is, however, quite 
clear from the neutron matter results (y p =0.0) that the 
isospin contents of the nucleons leads to higher radii dif- 
ferences at higher densities. Moreover, proton bags be- 
come larger than neutrons in medium, as matter goes 
denser. The neutron radius does not change much with 
isospin and for symmetric matter and neutron matter 
QMC and QMC<5 neutron radii almost overlap. 

In Fig. 2] we show the effective mass difference M* — 
M* for both QMC<5 and NL3<5, for y p =0.3 and 0.1. For 
NL3, we show the mass difference multiplied by a fac- 
tor of 0.1 in order to compare with QMC<5. The most 
striking result is the factor of ten difference between the 
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FIG. 3: (Color online) Radii ratios for neutron and proton, 
for symmetric matter and for y p —0.0 in the QMC5 model. 
Here we have taken R=0.6 fm~ 3 . 
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FIG. 4: (Color online) Effective mass differences M p — M n , for 
QMC<5 (blue, dotted) and NL3<5 (light green, dot-dashed), for 
different proton fractions: y p =0.1 (a) and y p =0.3 (b). Notice 
that the values for NL35 have been scaled by a 0.1 factor. 



p, n mass splitting in NL3 and QMC. In both models, 
this effect increases with baryon density, but it is worth 
remarking that the effective mass is larger for neutrons 
than for protons at lower densities, which is represented 
by the negative values in the figures. As referred before, 
this occurs because the proton and neutron masses were 
considered different at zero density. The effect of the 6- 
meson is to increase the M* — M* difference as it occurs 
in NL3<5 and other relativistic mean field models and in 
contrast to Brueckner-Hartree-Fock calculations [2l|, [22j 
In addition, for the same proton fraction, the crossing of 
the p, n effective mass curves (equal effective p, n masses) 
does not occur for the same densities as the proton and 
neutron radius, R p and R n (Fig. [3]). 



B. Instabilities 

In the present subsection we discuss the results for the 
instability region at subsaturation densities with QMC, 
QMC5 and the hadron models we have considered. In 
Fig. [5] we plot the spinodal curves for np matter. As 
referred before, they are defined by the points, for a given 
temperature, density and isospin asymmetry, that make 
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(Color online) Spinodal (thermodynamical instability 
for QMC, QMC<5, NL3, NL35 and TW. 



the curvature matrix of the free energy vanish. 

Both QMC and QMC5 present larger instability re- 
gions than NL3, NL35 and TW on the isoscalar direc- 
tion, p p — p n . This is possibly due to the a contribu- 
tion in the Lagrangian for QMC (both with and without 
5): the fields here attain magnitudes so as to minimize 
the bag energy, whereas in hadron models their mean- 
field values are determined by solving the relevant set of 
equations where nonlinearities show explicitly or through 
density dependent couplings. The extension of the spin- 
odal for p p = p n defines the density p s , corresponding 
to the density value for which the pressure of symmetric 
nuclear matter has a minimum and the incompressibility 
is zero [23]. We have included the values of this density 
for the different models in Table HI Recently [19( it was 
shown that within a Brueckner-Hartree-Fock calculation 
the isoscalar extension of the spinodal was much larger 
than NLWM and Skyrme interaction predictions. 

For large isospin asymmetries the presence of S reduces 
the instability region both in QMC and NL3. NL3 has 
a higher symmetry energy than NL35 for densities rang- 
ing from up to ~ po, and the same occurs for QMC 
although the differences are smaller. This means that 
highly asymmetric matter is less bound in NL3 and QMC 
than in NL35 and QMC<5. A mere inspection of the sym- 
metry energy is not enough to account for the differences 
in the instability region if different models are considered 
[23l ]. With the introduction of the <5-meson the scalar 
channel is not affected, and therefore we are essentially 
changing the isovector channel. We also notice that at 
Pp = Pn the curvature of the QMC spinodal is inter- 
mediate between NL3 and TW. As shown in [23| this 
curvature is defined by the symmetry energy and its first 
and second derivatives. 



The nuclear liquid-gas coexistence phase is character- 
ized by different isospin contents for each phase, i.e., the 
clusterized regions are more isospin symmetric than the 
surrounding nuclear gas, the so-called isospin distillation 
[ill |24| . The extension of the distillation effect is model 
dependent and it has been shown that NL3 and other 
NLWM parametrizations lead to larger distillation effects 
than the density dependent hadron models [g, |8|, [23j . On 
the other hand, the distillation effect was also studied 
within a BHF calculation [19( and a smaller distillation 
effect was generally obtained. 
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FIG. 6: (Color online) Direction of instability (eigenvector for 
negative eigenvalue A_) for j/ p =0.05 (a), 0.1 (b) and 0.3 (c). 
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In Fig. [6] we show the ratio of the proton versus the 
neutron density fluctuations corresponding to the unsta- 
ble mode. This ratio defines the direction of instability 
of the system. We show the results for different proton 
fractions (including rather small values), for the sake of 
studying the effectiveness of the models in restoring the 
symmetry in the liquid phase. 

We first compare the three models not including the 
5-meson: QMC, NL3 and TW. We see that QMC has a 
behaviour which is intermediate between NL3 and TW. 
The distillation effect for densities above 0.02 fm~ 3 is 
larger than the prediction of TW, but for the larger den- 
sities it also shows a tendency to decrease, contrary to 
NL3. The presence of 5-meson makes the distillation ef- 
fect more efficient. 

In Fig. [7] we plot the proton-neutron density fluctu- 
ation ratio as a function of the isospin asymmetry for 
a fixed nuclear density, p—0.06 fm~ 3 . We compare all 
models under study, NL3, NL3<5, TW, QMC and QMC<5, 
and include also the results derived from the BHF ap- 
proach [n| referred above. All the relativistic models 
predict larger fluctuation ratios than the corresponding 
value of Pp/pn, dotted line. The behaviour gives rise to a 
distillation effect, which, as referred before, is larger for 
NL3 and smaller for TW, with QMC presenting interme- 
diate values. The (5-meson stresses the distillation effect, 
clearly seen both in QMG5 and NL3(5. 

However, except for the very asymmetric matter (y p < 
0.02), both QMC and the other relativistic models pre- 
dict fluctuations with larger proton fractions than BHF. 
The instability properties within the BFH calculation of 
[l9| show also differences for the spinodal surface: the 
unstable region is larger, extending to larger densities 
both in the isoscalar an isovector directions and the cur- 
vature of the spinodal at p p — p n is much larger than 
the one of all the relativistic models considered. A larger 
extension of the unstable region is justified because the 
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FIG. 7: (Color online) Proton-neutron density fluctuation ra- 
tio versus the isospin asymmetry for a fixed nuclear density, 
p=0.06 fm" 3 . The BHF results were obtained by Vidana and 
Polls 0. 



saturation density is larger, 0.182 fm~ 3 . The shape of 
the spinodal itself depends on the density dependence of 
the symmetry energy and its derivatives. It would be im- 
portant to identify the properties that define the shape 
of the spinodal for the more asymmetric matter. 

We have studied subsaturation nuclear instabilities for 
both symmetric and asymmetric matte within the QMC 
model, with and without the inclusion of the <5-meson. In 
this model the nucleons are described as non-overlapping 
bags. We propose a parametrization for QMC<5 with a 
the symmetry energy slope value L = 102 MeV within 
the interval L = 88 ± 25 MeV proposed in p], 03 and 
which was determined from nuclear laboratory data. It 
was interesting to notice that for QMC and QMC<5 the 
quantity K asy defined in pd| . and which can be directly 
extracted from measurements of the isotopic dependence 
of the giant monopole resonances, falls inside the interval 
predicted by experiments. The BHF calculation predicts 
a value of K asy which is not very far, though lower than 
the interval obtained from GMR, K asy = 550±100 MeV. 

A comparison was done with the results obtained from 
a NLWM parametrization (NL3), one density dependent 
relativistic model (TW), along with BHF with the Ar- 
gonne V18 potential calculation. It was shown that the 
restoration of isospin symmetry, obtained by a distilla- 
tion effect, was more efficient in QMC with respect to 
TW but less efficient when compared with NL3. The 
spinodal surface within QMC is closer to TW although 
with a larger curvature at p p = p n and a slightly smaller 
instability extension at larger asymmetries, while the in- 
clusion of the <5-meson (QMC5) shrinks the asymmetric 
parts of the instability envelope. The BHF results of 
although with similar general properties, differ in the ex- 
tension and shape of the spinodal and in the amount of 
distillation predicted. 

A study of the QMC instability properties at finite 
temperatures is now under preparation. It is also im- 
portant to identify how the density dependence of the 
symmetry energy determines the shape of the spinodal 
for large isospin asymmetries. 
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